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Koide's mass formula is an empirical relation among the charged lepton masses which holds with 
a striking precision. We propose a mechanism for cancelling the QED correction to Koide's formula. 
This is discussed in an effective theory with (7(3) family gauge symmetry and a scenario in which 
this symmetry is unified with SU(2)l symmetry at f0 2 -10 3 TeV scale. 
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Among various properties of elementary particles, the 
spectra of the quarks and leptons exhibit unique pat- 
terns, and their origin still remains as a profound mys- 
tery. Koide's mass formula is an empirical relation among 
the charged lepton masses given by [l[ 
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which holds with a striking precision. In fact, substitut- 
ing the present experimental values of the charged lepton 
masses [2[ , the formula is valid within the present experi- 
mental accuracies. The relative experimental error of the 
left-hand side of eq. |T]) is of order 10~ 5 . 

Given the remarkable accuracy with which eq. ([T]) 
holds, there have been many speculations as to existence 
of some physical origin behind this mass formula [J, Q • 
Despite the attempts to find its origin, so far no realis- 
tic model or mechanism has been found which predicts 
Koide's mass formula within the required accuracy. The 
most serious problem one faces in finding a realistic model 
or mechanism is caused by the QED radiative correction. 
Even if one postulates some mechanism at a high energy 
scale that leads to this mass relation, the charged lepton 
masses receive the 1-loop QED radiative correction given 
by 
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where ffi(pi) and m pole denote the running mass defined in 
the modified - minimal-subtraction scheme (MS scheme) 
and the pole mass, respectively; [i represents the renor- 
malization scale. It is the pole mass that is measured 
in experiments. Suppose rhi(fx) (or the corresponding 
Yukawa couplings j/i(/i)) satisfy the relation ([T]) at scale 
U Mw Then m polc do not satisfy the same relation 
||: Eq. ([T]) is corrected by approximately 0.1%, which 
is 120 times larger than the present experimental error. 
Note that this correction originates only from the term 
— 3a/(47r) x fhi log(m?) of eq. |2]), since the other terms, 
which are of the form const, x fhi, do not affect the rela- 
tion |lj. We also note that log(mf) results from the fact 
that fhi plays a role of an infrared cut-off in the loop 
integral. 

The 1-loop weak correction is of the form const, x fhi 
in the leading order of fnf/M w expansion; the leading 



non-trivial correction is 0{Gpfh^ /ir) whose effect can 
be safely neglected. Other known radiative corrections 
are also negligible. 

Thus, if there is indeed a physical origin to Koide's 
mass formula at a high energy scale, we need to account 
for a correction to the relation {1} that cancels the QED 
correction. Since such a correction is absent up to the 
scale of 0(Mw) to our present knowledge, it must orig- 
inate from a higher scale. Then, there is a difficulty in 
explaining why the size of such a correction should coin- 
cide accurately with the size of the QED correction which 
arises from much lower scales. Up to date, no mechanism 
has been proposed to solve this problem. 

Among various existing models which attempt to ex- 
plain origins of Koide's mass formula, we find a class of 
models particularly attractive [6|. These are the models 
which predict the mass matrix of the charged leptons to 
be proportional to the square of the vacuum expectation 
value (VEV) of a scalar field (we denote it as $) written 
in a 3-by-3 matrix form: 
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Thus, (y/m e , y/nij,, yJm T ) is proportional to the diagonal 
elements of ($) in the basis where it is diagonal. The 
VEV ($) is determined by minimizing the potential of 
scalar fields in each model. Hence, the origin of Koide's 
formula is attributed to the specific form of the potential 
which realizes this relation in the vacuum configuration. 
Up to now, no model is complete with respect to sym- 
metry: Every model requires either absence or strong 
suppression of some of the terms in the potential (which 
are allowed by the symmetry of that model), without 
justification. 

In this paper, we propose a possible mechanism for 
cancelling the QED correction to Koide's mass formula, 
in the context of a theory with family (horizontal) gauge 
symmetry. In our study we adopt the mechanism eq. ([3]) 
for generating the charged lepton masses at tree level, due 
to the following reasons. First, models with this mecha- 
nism arc suited for pcrturbative analyses. Secondly, since 
$ is renormalized multiplicatively, the structure of radia- 
tive corrections is simple. 

We consider an effective theory with family gauge sym- 
metry which is valid up to some cut-off scale denoted by 
A (3> Mw)- Within such an effective theory, the charged 
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lepton masses may be induced by a higher-dimensional 
operator 



O 



A 2 



ipU $ik ®kj <P e Rj 
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Here, ipxA = { v Li:GLi) T denotes the left-handed lepton 
SU(2)l doublet of the z-th generation; e R j denotes the 
right-handed charged lepton of the j-th generation; if de- 
notes the Higgs doublet field; $ is a 9-component scalar 
field and is singlet under the Standard Model (SM) gauge 
group. We suppressed all the indices except for the gen- 
eration (family) indices i,j,k= 1, 2, 3. (Summation over 
repeated indices is understood throughout the paper.) 
The dimensionless Wilson coefficient of this operator is 
denoted as Once $ acquires a VEV, the operator 

O will effectively be rendered to the Yukawa interactions 
of the SM; after the Higgs field also acquires a VEV, 
(<p) = (0,v cw /V2) T with v cw ~ 250 GeV, the operator 
will induce the charged-lepton mass matrix of the form 
eq. at tree level: 



mY cc = 



V2A 2 
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We consider radiative corrections to the above mass 
matrix by the family gauge interaction. First we con- 
sider the case, in which the gauge group is SU(3) and 
both ipL and e R are assigned to the 3 (fundamental rep- 
resentation) of this symmetry group. With this choice 
of representation, however, Koide's formula is subject to 
a severe radiative correction unless the family gauge in- 
teraction is strongly suppressed. In fact, the 1-loop cor- 
rection by the family gauge bosons induces an effective 
operator 



V ~ — x k ipLi <p em x 
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hence corrections universal to all the charged-lepton 
masses, (5m e , Sm^, 5m T ) oc (1, 1, 1), are induced. This is 
due to the fact that the dimension-4 operator ipn ipem 
is not prohibited by symmetry. Here, otp — <?^/(47r) de- 
notes the gauge coupling constant of the family gauge in- 
teraction. As noted above, corrections which are propor- 
tional to individual masses do not affect Koide's formula; 
oppositely, the universal correction violates Koide's for- 
mula strongly. In order that the correction to Koide's for- 
mula cancel the QED correction, a naive estimate shows 
that qf/tt should be order 10~ 5 , provided that the cut- 
off A is not too large and that the above operator O' is 
absent at tree level. If O' exists at tree level, there should 
be a fine tuning between the tree and 1-loop contribu- 
tions. The situation is similar if the family symmetry is 
0(3) and both ipL and e R are in the 3. In these cases 
we were unable to find any sensible reasoning for the 
cancellation between the QED correction and the correc- 
tion induced by family gauge interaction, other than to 
regard the cancellation as just a pure coincidence. Hence, 
we will not investigate these choices of representation fur- 
ther. 



In the case that ipL is assigned to 3 and e R to 3 (or vice 
versa) of U{3) family gauge group, (i) the dimension-4 
operator tpLifem is prohibited by symmetry, and hence 
corrections universal to all the three masses do not ap- 
pear; and (ii) marked resemblance of the radiative cor- 
rection to the QED correction follows. We show these 
points explicitly in a specific setup. 

We denote the generators for the fundamental repre- 
sentation of 17(3) by T a (0 < a < 8), which satisfy 



tr (T a T ) = i 6 a ^ 



rj-iCt fjl t 
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T° = ^=1 is the generator of U(l), while T a (1 < a < 8) 

are the generators of SU(3). 

We assign ipL to the representation (3,1), where 3 
stands for the SU(3) representation and 1 for the U(l) 
charge, while e R is assigned to (3, —1). Under U(3), the 
9-component field $ transforms as three (3, l)'s. Explic- 
itly the transformations of these fields are given by 



i'L^U^L, e R ^U*e R , $ -»• U 3> 



(8) 



with U = exp (i9 a T a ). 

We assume that the charged-lepton mass matrix is in- 
duced by a higher-dimensional operator O^' similar to 
O in eq. Q. We further assume that (5>) can be brought 
to a diagonal form in an appropriate basis. Thus, in this 
basis after $ and <p acquire VEVs, turns to the 

lepton mass terms as 
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■ ipL ^d(^) 2 e R , 



(9) 



where 



vi{fi) 
v»(ji) 
v 3 ( P ) 



Vi{n)>0. (10) 



When all Vi are different, J7(3) symmetry is completely 
broken by ($) = &d, and the spectrum of the U(3) gauge 
bosons is determined by <&d. 

Note that the operator O in eq. is not invariant 
under the U(3) transformations eq. ([8]). As an example 
of O", one may consider 



0^ = ^fc***W 



(11) 



It is invariant under a larger symmetry U{3) x 0(3), un- 
der which <I> transforms as $ — * U$0 T (0 T = 1). In 
this case, we need to assume, for instance, that the 0(3) 
symmetry is gauged and spontaneously broken at a high 
energy scale before the breakdown of the U{3) symme- 
try, in order to eliminate Nambu-Goldstone bosons and 
to suppress mixing of the J7(3) and 0(3) gauge bosons. 
In any case, the properties of 0^ given by eqs. ([9]) and 
pop are sufficient for computing the radiative correction 
by the U(3) gauge bosons to the mass matrix, without 
an explicit form of . 
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We take the U(l) and £{7(3) gauge coupling constants 
to be the same: 



a U(l) = a SU(3) — a F ■ 



(12) 



We compute the radiative correction in Landau gauge, 
which is known to be convenient for computations in 
theories with spontaneous symmetry breaking. Through 
standard computation, one obtains 



5m 



pole 
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log 



Vi(n) 2 . 



m t (/x) , (13) 
(14) 



Here, c is a constant independent of i. The Wilson co- 
efficient n^(fi) is defined in MS scheme. Vi(fj.) are de- 
fined as follows: The VEV of <E> at renormlization scale fj,, 
&d{[J-) = (^0-0) given by eq. (|10p . is determined by min- 
imizing the 1-loop effective potential in Landau gauge 
(although we do not discuss the explicit form of the ef- 
fective potential); $ is renormalized in MS scheme. We 
ignored terms suppressed by m 2 /v 2 {<^L 1) in the above 
expression. Note that the pole mass is renormalization- 
group invariant and gauge independent. Therefore, the 
above expression is rendered gauge independent if we ex- 
press Vi(fj) in terms of gauge independent parameters, 
such as coupling constants defined in on-shell scheme. 

The form of the radiative correction given by eqs. (I13p 
and (|14[) is constrained by symmetries and their breaking 
patterns. As the diagonal elements of the VEV, v$ > 
v-2 > vi > 0, are successively turned on, gauge symmetry 
is broken according to the pattern: 



U(3) -> 17(2) -> £7(1) -> nothing. 



(15) 



At each stage, the gauge bosons corresponding to the 
broken generators acquire masses and decouple. Further- 
more, the vacuum <i>d and the family gauge interaction 
respect a global U(l) 3 symmetry generated by 



^Pl^V^ l , e R -> V* e R , $ d -> V $ d V* , (16) 



with 



V = 










e i<i>2 





e r. 



(17) 



Although after symmetry breakdown, eq. ©, is not 
invariant under this transformation, the variation can be 
absorbed into a redefinition of u,. As a result, the lepton 
mass matrix has a following transformation property: 



M, 



VM e V* 



(18) 



Vi — *Vi cxp(i0i) 



This is satisfied including the 1-loop radiative correction. 
The symmetry breaking pattern eq. (|15[) and the above 



transformation property constrain the form of the radia- 
tive correction to <5m? ole oc v 2 [log(|wi| 2 ) + const.], where 
the constant is independent of i. Note that \vi\ 2 in the 
argument of logarithm originates from the gauge boson 
masses, which are invariant under «j — > Vi exp(i^). 

The universality of the SU(3) and 17(1) gauge cou- 
plings eq. (|12[) is necessary to guarantee the above sym- 
metry breaking pattern eq. (| 15[) . One may worry about 
validity of the assumption for the universality, since the 
two couplings are renormalized differently in general. 
The universality can be ensured approximately if these 
two symmetry groups are embedded into a simple group 
down to a scale close to the relevant scale. There are more 
than one ways to achieve this. A simplest way would be 
to embed SU(3) x U(l) into £{7(4). It is easy to ver- 
ify that the 4 of SU(4) decomposes into (3, --|) (1, |) 
under SU(3) x U(l). Hence, the 6 (second-rank anti- 
symmetric representation) and 6 of SU(4), respectively, 
include (3, -1) and (3, 1). 

Within the effective theory under consideration, the 
QED correction to the pole mass is given just as in eq. (|2|) 
with fhi(ii) replaced by rrii(p). Recall that corrections 
of the form const, x m, do not affect Koide's formula. 
Noting \ogvf = ifogm 2 + const., one observes that if 
a relation between the QED and family gauge coupling- 
constants 
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(19) 



is satisfied, the 1-loop radiative correction induced by 
family gauge interaction cancels the 1-loop QED correc- 
tion to Koide's mass formula. 

Suppose the relation (TT9")) is satisfied. Then 



mf ole oc Vi{nY 



(20) 



holds with a good accuracy. This is valid for any value 
of /i. This means, if Vi(jj) satisfy 



vMm) 2 +v 2 (fi) 2 +v 3 (/i) 2 



(21) 



at some scale /x, Koide's formula is satisfied at any scale 
/j. This is a consequence of the fact that $ is multiplica- 
tively renormalized. Generally, the form of the effective 
potential varies with scale [i. If the relation (|2T|) is real- 
ized at some scale as a consequence of a specific nature 
of the effective potential (in Landau gauge) , the same re- 
lation holds automatically at any scale. Although these 
statements are formally true, physically one should con- 
sider scales only above the family gauge boson masses, 
since decoupling of the gauge bosons is not encoded in 
MS scheme. An interesting possibility is to use eq. (|20|) 
to relate the charged lepton pole masses with the VEV 
at the cut-off scale, i.e. /x = A, which sets a boundary 
(initial) condition of the effective theory. 

The advantages of choosing Landau gauge in our com- 
putation are two folds: (1) The computation of the 1- 
loop effective potential becomes particularly simple (as 
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well known in computations of the effective potential in 
various models); in particular there is no 0(otp) cor- 
rection to the effective potential. (2) The lepton wave- 
function renormalization is finite; as a consequence, the 
Pfj.1 11 part of the lepton self-energy is independent of gen- 
eration. Due to the former property, there is no 0(o.f) 
correction to the relation eq. (f2"Tj) if it is satisfied at 
tree level. Due to the latter property, the correction to 
Koide's formula is determined by renormalization of O™' 
alone, and a simple relation to follows. 

Let us comment on gauge dependence of our predic- 
tion. If we take another gauge and express the radiative 
correction Sm^ ole in terms of (<£>(/j)), the coefficient of 
log(/i 2 /($) 2 ) changes, and other non-trivial flavor de- 
pendent corrections are induced. Suppose the relation 
eq. ([21]) is satisfied at tree level [nj]. The VEV ($) in an- 
other gauge receives an O(ttf) correction, which induces 
a correction to eq. (|21[) at 0(ap)- These additional cor- 
rections to <5m? olc at 0(cxf) should cancel altogether if 
they are reexpressed in terms of the tree-level v^s which 
satisfy eq. (|2"Tj) . since the 0(ap) correction to the rela- 
tion (|2ip vanishes in Landau gauge. General analyses on 
gauge dependence of the effective potential may be found 
in Q. 

Now we speculate on a possible scenario how the rela- 
tion (p~9]) may be satisfied. The scale of a is determined 
by the charged lepton masses, while the scale of cxf is de- 
termined by the family gauge boson masses, which should 
be much higher than the electroweak scale. Since the rel- 
evant scales of the two couplings are very different, we 
are unable to avoid assuming some accidental factor (or 
parameter tuning) to achieve this condition. Instead we 
seek for an indirect evidence which indicates such an ac- 
cident has occurred in Nature. The relation (fT9|) shows 
that the value of qj? is close to that of the weak gauge 
coupling constant aw, since sin 2 9w{Mw) is close to 1/4. 
In fact, within the SM, j otw{n) approximates a(m T ) at 
scale [i ~ 10 2 -10 3 TeV. Hence, if the electroweak SU{2) L 
gauge group and the U(3) family gauge group are uni- 



fied around this scale, naively we expect that a ~ \ocf 
is satisfied. Since aw runs relatively slowly in the SM, 
even if the unification scale is varied within a factor of 3, 
Koide's mass formula is satisfied within the present ex- 
perimental accuracy. This shows the level of parameter 
tuning required in this scenario. 

We assume that anomalies introduced by the cou- 
plings of fermions to family gauge bosons are cancelled, 
which requires existence of fermions other than the SM 
fermions. Furthermore, we assume that all the additional 
fermions acquire masses of order ($) or larger after the 
spontaneous breakdown of U(3), so that they decouple 
from the SM sector at and below the electroweak scale. 

A characteristic prediction of the present scenario is 
the existence of lepton flavor violating processes at 10 2 - 
10 3 TeV scale. For instance, assuming that the down- 
type quarks are in the same representation of U(3) as 
the charged leptons, and that the mass matrices of 
the charged leptons and down-type quarks are simul- 
taneously diagonalized in an appropriate basis, we find 
T(Kl — > e/i) m 2 r7i^/|./(167rw|). Comparing to the 
present experimental bound Bt(Kl — > e/i) < 4.7 x 10~ 12 
[2fl , we obtain a limit v 2 >5x 10 2 TeV. Naively this limit 
may be marginally in conflict with the estimated unifica- 
tion scale in the above scenario. This depends, however, 
rather heavily on our assumptions on the quark sector. In 
the case that there exist additional factors in the quark 
sector which suppress the decay width by a few orders 
of magnitude, we may expect a signal for Kl — > efj, not 
far beyond the present experimental reach. Predictions 
concerning purely leptonic processes are less model de- 
pendent, but expected event rates are far below present 
experimental sensitivities. 

Models, which predict a realistic charged lepton spec- 
trum incorporating the mechanism proposed in this pa- 
per, will be discussed elsewhere floT ]. 
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